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Abstract 

We introduce a geometric generalization of the 0(A)-field theory that describes 
A^-colored membranes with arbitrary dimension D. As the 0(A)-model reduces 
in the limit N — > to self-avoiding polymers, the A-colored manifold model leads 
to self-avoiding tethered membranes. In the other limit, for inner dimension D — ► 
1, the manifold model reduces to the 0(A)-field theory. We analyze the scaling 
properties of the model at criticality by a one-loop perturbative renormalization 
group analysis around an upper critical line. The freedom to optimize with respect 
to the expansion point on this line allows us to obtain the exponent v of standard 
field theory to much better precision that the usual 1-loop calculations. Some other 
field theoretical techniques, such as the large N limit and Hartree approximation, 
can also be applied to this model. By comparison of low and high temperature 
expansions, we arrive at a conjecture for the nature of droplets dominating the 
3cf-Ising model at criticality, which is satisfied by our numerical results. We can 
also construct an appropriate generalization that describes cubic anisotropy, by 
adding an interaction between manifolds of the same color. The two parameter 
space includes a variety of new phases and fixed points, some with Ising criticality, 
enabling us to extract a remarkably precise value of 0.6315 for the exponent v in 
d = 3. A particular limit of the model with cubic anisotropy corresponds to the 
random bond Ising problem; unlike the field theory formulation, we find a fixed 
point describing this system at 1-loop order. 
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1 Introduction 



Field theoretical models are particularly suited for description of universal quantities 
which do not depend on the details of the system. This is epitomized by systems un- 
dergoing symmetry breaking continuous phase transitions. Their critical behavior is de- 
scribed by a set of exponents which are completely characterized by dimension and the 
underlying symmetry (the number of components of the order parameter). Universality 
is assured since the microscopic details are averaged out, and do not effect the large scale 
fluctuations (for a review, see Ref. [1]). The universal features are thus captured by the 
O(N) model, which is a field theory for the statistics of iV-component spins. A variety 
of techniques have been developed to obtain the critical behavior of this model; possibly 
the most successful is the renormalization group procedure [1] which analytically justifies 
the concept of universality. The most technically convenient implementations are field 
theoretical methods, e.g. the ^-expansion about the upper critical dimension of 4, an 
expansion about the lower critical dimension of 2, and exact resummations in the large 
N limit. (For an overview of these techniques, see Ref. [2].) The best studied method 
is the e-expansion about the upper critical dimension 4, where calculations have been 
performed up to 5-loop order. Together with resummation techniques which take care of 
the large-order behavior known from instanton calculus, this is a very powerful tool for 
extracting critical exponents. 

On the other hand, field theories have strong connections to geometrical problems 
involving fluctuating lines. For example, the motion of particles in space-time describes a 
world-line. Summing over all world-lines, weighted by an appropriate action, is the Feyn- 
man path integral approach to calculating transition probabilities, which can alternatively 
be obtained from a quantum field theory. Another example is the high-temperature ex- 
pansion of the Ising model. The energy-energy correlation function can be expressed 
as a sum over all self-avoiding closed loops which pass through two given points. The 
generalization to iV-component spins is straightforward: The partition function of the 
corresponding 'loop model' is obtained by summing over all configurations of a gas of 
closed loops, where each loop comes in N colors, or has a fugacity of N. In the limit 
N — > 0, only a single loop contributes, giving the partition function of a closed self- 
avoiding polymer. For N > the model describes polymers which can break up and 
polymerize dynamically like liquid sulfur [3-5]. 

A more direct approach to study self-avoiding polymers was developed by Edwards 
and Des Cloizeaux [6-8] . In this approach, hard self- avoidance is replaced by a soft short 
range repulsive interaction between the monomers. The repulsive interaction is then stud- 
ied perturbatively by expanding about ideal random walks. Here too, the perturbative 
expansion can be reorganized into a renormalization group about the upper critical dimen- 
sion 4, which was shown [9] to be equivalent to the perturbation expansion of <^ 4 -theory 
in the limit N — > 0. This equivalence holds both on the formal and on the perturbative 
levels, providing two apparently different approaches for calculating the same exponents. 

There is much work in the field theory community on generalizing results for fluctu- 
ating lines to entities of other internal dimensions D. The most prominent example is 
the work on string theories, which describe D = 2 world sheets. An earlier example is 
provided by the correspondence between gauge theories and random surfaces [10,11]. The 



3 



low temperature expansion of the Ising model in <i-dimension also results in a sum over 
surfaces that are d — 1 dimensional. For d — 3, the surfaces are made out of plaquettes, 
the basic objects of lattice gauge theories. 

The simplest generalization of linear polymers is to "tethered" (or polymerized) sur- 
faces [12,13], which have a fixed internal connectivity, and are thus simpler than their 
gauge theory counterparts. For theoretical analysis, it is convenient to further general- 
ize to membranes of arbitrary (inner) dimension D, interpolating between polymers for 
D — 1 and membranes for D = 2. Simple power counting indicates that the self-avoiding 
interaction is relevant only for dimensions d < d c = AD/ '(2 — D), making possible an 
e = 2D — d{2 — D)/2 ~ (d c (D) — <f)-expansion, which was first carried to 1-loop order 
around this line in Refs. [14-18]. To obtain results for polymers or membranes, one now 
has the freedom to expand about any inner dimension D, and the corresponding upper 
critical dimension of the embedding space [19]. This freedom can be used to optimize the 
calculation of critical exponents. 

Following more rigorous analysis of this novel perturbation series [20-24], recently 

2- loop calculations were performed for membranes with inner dimension D between 1 
and 2 [25,26]. The results have been applied to polymers, where the swelling exponent v 
has been found to be 0.59 in an appropriate extrapolation scheme at both 1- and 2-loop 
order. By contrast, in standard field theory, the 1-loop result is an underestimate, while 
the 2-loop result is an overestimate by a similar amount. For self- avoiding membranes in 

3- dimensional space, 2-loop calculations predict an isotropic fractal phase with dimension 
of about 2.4 [25,26]. In addition, there have been extensive numerical studies [27-36], and 
a few experiments on graphite oxide layers [37-39]. 

In this article, we reverse the analogy that leads from the O(N) model to self- avoiding 
polymers: The idea is to generalize the high temperature expansion of the O(N) model 
from a gas of self-avoiding loops of fugacity N, to a similar gas of closed fluctuating 
manifolds of internal dimension D. The primary goal is to obtain a novel analytical 
handle on the field theory for D — 1, and we do not insist that the models for general D 
correspond to any physical problem. Given this caveat, the generalization is not unique. 
Encouraged by its success in polymer theory, we study the generalization to tethered 
manifolds, and in addition restrict ourselves to the genus of hyper-spheres. For this class 
of surfaces calculations are simpler; in particular yielding excellent values of the exponent 
v in the limit of polymers [25,26]. We have chosen hyperspheres as they have no additional 
anomalous correction exactly at D = 2. The resulting manifold theory depends on two 
parameters N and D, whose limiting behaviors reduce to well known models, as indicated 
in the diagram below. 

The rest of this paper is organized as follows: In Sec. 2 we review the high-temperature 
expansion for the O(N) model to explicitly indicate how it relates the sum over self- 
avoiding loops. The different renormalization group schemes used in the literature to 
study polymers may potentially lead to some confusion regarding the connection to y? 4 - 
theory. To clarify the situation, we explicitly compare and contrast the two main schemes 
in Sec. 3. We then generalize the model in Sec. 4 to polymers of N different colors, 
making the connection to the O(N) ^/-theory. Our central result, obtained in Sec. 5, is 
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Figure 1.1: Schematic description of the new model, and its limits. 



the first-order expansion for the exponent v(D, N, d), given by 
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The above result depends on a function c(D), which at the microscopic level is related to 
the relative strengths of self-avoidance between parts of the same manifold, and between 
different manifolds. In Sec. 6 we propose two choices, c(D) = 1, and c(D) = D, for this 
parameter. The numerical values resulting from several extrapolations are discussed in 
Sec. 7. For the O(N) model in d — 3, our best extrapolations for the exponent v are 
equal to 0.601, 0.646, 0.676, and 0.697, for JV = 0, 1, 2, and 3 respectively. 

The ambiguity associated with c(D) disappears in the N — > oo limit, where it is 
possible to exactly sum the dominant diagrams in the perturbation series. This result, 
along with a mean field variational estimate, is presented in Sec. 8. 

The low temperature expansion of the ci-dimensional Ising model (N = 1) provides 
another route to sums over random surfaces. As discussed in Sec. 9, the sum is over 
surfaces of internal dimension D — d— 1, embedded in d dimensions. However, in crucial 
difference with tethered manifolds, it is necessary to sum over all internal metrics (con- 
nectivities). For A^ — > 0, this difference is known to lead to quite drastic geometries. In 
particular, sums over a single random surface are dominated by singular configurations 
which in fact resemble branched polymers [40]. However, for N ^ 0,it may be entropically 
advantageous to break up a singular spike into many bubbles. If so, a description in terms 
of fluctuating hyper-spherical surfaces may not be too off the mark. The requirement that 
the dual high and low-temperature expansions of the Ising model partition function have 
the same singularity, leads to a putative identity a(N = l,d,D) = a(N = l,d,d— D). 
The numerical tests based upon the 1-loop result of Eq. (1.1) appear to support this 
conjecture. 

While obtaining better exponents for the 0(N) model is an important goal, our gen- 
eralized approach is more valuable if it also applies to other local field theories. The 
simplest extension of the A^-component y? 4 theory includes cubic anisotropy by adding a 
term proportional to J2 a ft- The geometrical interpretation of this term is an additional 
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interaction that only operates between manifolds of the same color. In Sec. 10 we develop 
the corresponding manifold extension, whose renormalization equations involve two inter- 
action parameters, leading to a variety of new phases and fixed points, some with Ising 
criticality. A particular scheme in this model yields an Ising exponent of v* = 0.6315 in 
d — 3, which is indistinguishable from the many loop calculations [2]! The consequences 
of this generalization for the random bond Ising model are explored in Sec. 11. In the 
standard description with D — 1 there is no fixed point for the random bond Ising model 
at 1-loop order, necessitating a \fe expansion. By contrast, we do find a fixed point at 
this order for D ^ 1. 

A number of technical discussions are relegated to the appendices: Appendices A and 
B present derivations of some properties of the renormalization group factors used in the 
text. For completeness, and convenience of the reader, some technical details of dealing 
with divergences in the perturbation series for D ^ 1 are presented in Appendices C and 
D. Finally, Appendix E deals with the question of what happens if manifolds of other 
topology are also considered. 



2 The 0(iV)-model in the high-temperature expansion 

In this section, we briefly review the high-temperature expansion of the O(N) model. 
(For more extensive reviews, see Refs. [11] and [41].) The Hamiltonian is 

H = - JN ]T Si ■ Sj , (2.1) 
where the sum runs over all nearest neighbors of a ci-dimensional cubic lattice. To obtain 

— # 

the partition function, we have to integrate over all Si subject to the constraint that 
|<Si| = 1, resulting in (K — f3J) 

e-W = ( TT e"*^ . (2.2) 
{Si} J{s t } J* 

The high-temperature expansion is obtained by expanding the exponential factors in Eq. 
(2.2) as 

e NKS t -s 3 = 1 + NK §, + . (2.3) 

Typically, only the first two terms in the Taylor-expansion are retained. This is justified 
as we are only interested in universal quantities, for which the weight is already not unique 
and may be modified [exp(NKSi ■ Sj) — ► 1 + NKSi ■ Sj] in order to cancel subsequent 
terms in the Taylor expansion. 

We can represent the various terms in the perturbation expansion in the following 
manner (see Refs. [11,42]): For each term NKSi ■ Sj, we draw a line connecting sites 
i and j. At any given site i, up to 2d such lines may join. The integral over the spin 

— * 

Si is non-zero, if and only if an even number of bonds end at site i. For calculational 
convenience, we normalize the integrals by the corresponding solid angle such that 



/ 



dSi = 1 . (2.4) 
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Figure 2.1: Some terms in the high-temperature expansion of the (9(iV)-model. 



Let us now study the first few terms in the perturbation expansion (see figure 2.1). The 
diagram (a) is 

(a) = (KN) 4 J ' dSi . . . dS 4 S%S%S%SigZSlS{S{ . (2.5) 
To do the integrations, note that 

J dSj 2 = J dSU = l, (2.6) 

and therefore 

dS, S?S? = . (2.7) 



N 

Performing all but the last integration in Eq. (2.5), we obtain 

(a) = K A N J dSi Si = K 4 N . (2.8) 
For any non-intersecting loop, this result is easily generalized to 

j^numbcr of links jy (2 9) 

i.e. every closed loop contributes a factor of N. Let us now analyze what happens when 
loops intersect and to this aim calculate configuration (b). Doing all but the integration 

— * 

over Si, we obtain 

(b) = K 8 N 2 J dSi (Sf) 2 = K 8 N 2 = (a) 2 . (2. 10) 

Two configurations which have one common point thus give the same contribution in the 
high-temperature expansion as if they were disjoint. This is not the case if they have 
one bond in common, see (c). The integral contains an odd power of the field Si, and 
therefore 

(c) = 0. (2.11) 

This high-temperature series can thus be reinterpreted as the sum over all self-avoiding 
(non- intersecting) loops. Bonds are totally self-avoiding, see e.g. configuration (c), while 
vertices are also partially self-avoiding as can be seen from the following argument. There 
are 3 possible ways to build up configuration (b): One may take 2 small loops, but there are 
also 2 possibilities to use one loop only. The latter have to be excluded from the partition 
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function. (There are additional constraints associated with multiple intersections.) On 
the other hand, as we are only interested in universal quantities, taking precise account 
of these configurations should be irrelevant as long as bond-self-avoidance is present. In 
the direct polymer approach of Edwards and Des Cloizeaux [6-8] discussed below, this 
corresponds to taking a smaller initial (bare) coupling constant. 

A single loop can now be viewed as a random walk, i.e. as the trace of a particle 
moving under Brownian motion. The corresponding Hamiltonian is 

H = [ L dxhvr(x)) 2 + Lt , (2.12) 
Jo 4 

where r(x) G JR d is the trajectory of the particle at time x (equivalently, x is the polymer 
arc-length). The total length of the loop is L — J dx. In addition, one has to demand 
that the particle returns to its starting-point, i.e. that the polymer is closed. To make 
it self-avoiding, Edwards and Des Cloizeaux [6-8] added an explicit repulsive interaction 
upon contact, leading to 

H= [ L dxh\/r(x)) 2 + ^- [ L dx [ L dy5 d (r(x) - r(y)) + Lt . (2.13) 
Jo 4 4 Jo jo 

The factors of 1/4, as well as the normalization hidden in 5, are chosen for convenience 
and will be explained later on; fi sets the renormalization scale. In the high-temperature 
expansion, there appear loops of all sizes. We thus have to sum over all different lengths 
of the polymer, weighted by a chemical potential t conjugate to the length, mimicking the 
constant K in Eq. (2.2). To avoid possible confusion, let us stress that although closely 
related, \nK and t are not identical. While K is defined as the fugacity for the length 
of the lattice walk, the chemical potential t is conjugate to the coarse-grained length. 
In principle, the same lattice walk can be represented by curves r(x) of different length 
L. However, as far as universal quantities are concerned, this is unimportant. Both 
parameters have to be tuned to reach the critical point, and only their deviations from 
the critical value, but not the critical value itself, have some physical correspondence. 



3 Renormalization group for polymers 

We now discuss the perturbation expansion of the Hamiltonian in Eq. (2.13). Let us 
start with the correlation functions of the free (non self-avoiding) polymer. One has to 
be careful in distinguishing between open and closed polymers which will be denoted by 
subscripts and c respectively. For open (or closed, but infinitely long) polymers, the 
correlation function 

C o (x) = ±(±(r(x)-r(0))^ o , (3.1) 
is the solution of the Laplace equation 

\aC (x) = 5(x) , (3.2) 

which is easily found to be 

C {x) = \x\ . (3.3) 
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For closed polymers, Eq. (3.3) has to be modified. The reason is that the information has 
two equivalent ways to travel around a polymer loop of size L, leading to 

C c (x) = N(L ~ N) , V |x| < L . (3.4) 

We next calculate the weight of a polymer of length L. For open polymers this is 
simply 

(3.5) 

where t is the chemical potential. For closed polymers, an additional factor of 

5 d (r(L) -r(0))) (3.6) 



has to be added, which measures the probability to find a closed polymer among all open 
polymers. The expectation value therefore is taken with respect to the weight for an open 
polymer, and calculated as follows: 

(5 d (r(L) - r(0))) o = ^ ^Hr(L)-r(o))^ 

-k 2 C (L) 



= /e 

Jk 



lk 



e -k*L 



k 

= L~ d ' 2 . (3.7) 
The normalizations of 5 d and f k are chosen for calculational convenience such that 



lk 



The same normalizations are also used to incorporate self-avoidance as discussed later. 

To get the quantities obtained in the high-temperature expansion of the loop model 
introduced above, we still have to integrate over all possible lengths of the polymer. We 
define the partition function for a single polymer as 



(3.9) 



We have chosen to integrate over a logarithmic scale (JdL/L) in order to make the 
integration measure dimensionless. The factor L counts for the number of points which 
may be taken as origin. Our final check, however, is that we obtain the same result as in 
the free (Gaussian) field theory. 

Additional insight is obtained from a different way to calculate zf\ If we do not 
perform the last integral in Eq. (3.7), Eq. (3.9) becomes 

Z<°»= /die-" j e-^ 
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As suggested graphically above, this term of the polymer-perturbation theory is equivalent 
to a term in the perturbation theory of the field-theoretical description of the 0(iV)-model. 
In the usual treatment of the O(N) field theory, the hard constraint of \S\ — 1 is replaced 
in favor of a soft constraint, implemented by the Hamiltonian 



n 



O(N) 



d d r 



\{VS{r)Y + \s\r) + h J^{S\r)f 



(3.11) 



In this description, one has to take the limit iV — > in order to allow for only 1 connected 
piece. (Remember that every closed loop counts a factor of N.) This equivalence, first 
pointed out by De Gennes [9], is not accidental and can be proven both perturbatively 
and by formal manipulations of the functional integral [2]. It reflects the fact that both 
the field-theoretic formulation of the 0(iV)-model, as well as its lattice equivalent, belong 
to the same universality class. The reader not familiar with this statement is invited to 
manipulate a few other terms in the perturbation expansion. In the following discussion, 
we shall demonstrate this equivalence for all diagrams encountered. 

We now perform the perturbation expansion of the polymer Hamiltonian in Eq. (2.13). 
The first term is the expectation value of one ^-interaction with respect to the free theory 
of a closed polymer, integrated over all positions of the interaction on the polymer of 
length L, and then over all polymer-lengths. This is explicitly 



poo , pL pL p 

/ dLL- d / 2 e~ Lt / dx dy 

JO JO Jo Jk 



dLe 



-Lt 



dx / dy 



\x 



-r(v))\ 

I c 

— \x 



v\) 



L 

p OO P J-j PX 

= 2 / dLe~ Lt / dx dz [z(L - z)]- d/2 
Jo Jo Jo 

I f 

dy > e -t(x'+y'+z) 

Jp 



-d/2 



L -d/2 



dz 



dx' 



p 2 z 



Ik 



-k 2 (x'+y') 



= 2 If 

Jk Jp 



p (k 2 +tfp 2 + t 

o o 



(3.12) 



The relation to (^ 4 -theory is again apparent: The integrals in Eq. (3.12) are ultra-violet 
divergent. The leading divergence is subtracted via a finite part prescription, the sub- 
leading term is treated via dimensional regularization as a pole in 



e = 2 - d/2 . 



(3.13) 



(Note the factor of 2 difference from the more usual definition of e — 4 — d.) 

Let us now introduce a renormalized Hamiltonian. Three renormalizations may be 
required: A renormalization of the field r, of the coupling constant b, and of the chemical 
potential t. Denoting the bare quantities with a subscript o, we set 



r = \[Zr , 

b = jj £ Z d / 2 Z b b . 



(3.14) 
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This yields the renormalized Hamiltonian 

H = Z jdx^(Vr(x)) 2 + b -^ jdx jdy~5 d (r(x)-r(y)) + Z t t Jdx , (3.15) 

where fx sets the renormalization scale. It is possible to subtract at the scale of the 
renormalized chemical potential t, but this turns out to be rather confusing when deriving 
the renormalization group equations. We can now eliminate the divergence in Eq. (3.12) 
by setting 

Z t = l- ^Res(Q) . (3.16) 

This is seen by expanding e~ n with H given in Eq. (3.15). From Eqs. (3.9) and (3.10), 
we read off the numerical value of , yielding 

Z t = l + ^. (3.17) 

The next step is to study the renormalization of the interaction, to which the following 
two diagrams contribute 

, and j V{ . (3.18) 



X 



To calculate the first diagram, change coordinates to x and y , which indicate the points 
midway between the contacts on each polymer. The shorter relative distance between 
these points on each polymer is denoted by x (or y), while the longer one is indicated 
by fl x (or fly). The arbitrariness in this choice leads to a combinatorial factor of 2 per 
polymer loop, for an overall coefficient of 4. For each contribution of 



X 



f f e -t(n x +n y +x+y) / e i(f£+fc)(r(*o + f )-r(yo+%)) e i(^-k)(r(x -%)- r (y -%)) 

Jk ,x ,y Jk,x,y \ 

(3.19) 



short distance singularities appear in the integration over x and y. The leading term in 
the short distance expansion is 

f /ik (r(x )-r(y ))\ e -t(n x +n y ) f Q -k 2 (C c (x)+C c (y)) e -t(x+y) _ (3.20) 

Jk ,x ,yo Jk,x,y 

For small arguments, the correlation function can be approximated by its infinite volume 
limit, leading up to sub leading terms to 



The final result is 



X 



+ subleading terms . (3.22) 
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The second diagram in Eq. (3.18) has already appeared in Eq. (3.12), and we can sym- 
bolically write 




(DO O 



k (k 2 + t) 2 



r U - ^ t 2 - d ' 2 



(3.23) 



This diagram appears with a combinatorial factor of 2 for its left-right asymmetry, and 
another factor of 2 for the possibilities to put the single point on i \ j . 

Adding these contributions yields the following renormalization factor at 1-loop order 
(note that the combinatorial factors of 4 cancel with that of the b/A in the Hamiltonian 
Eq. (3.15)), 

». iy^^ r] ,.-.>/■ (3 24) 

No field renormalization is necessary (Z = 1). 

The next step is to calculate the renormalization group functions, which measure the 
dependence of the renormalized quantities upon a change of the renormalization scale /x, 
while keeping the bare values fixed. The derivation of these functions is given in appendix 
A, and results in a so-called /3-function 



z b = l + -Res ( 



^=4 



and a scaling function for the field R 



-eb 



l + &£lnZ 6 + f&£lnZ ' 



(3.25) 



^) = i->)|ln(^ t ). 



(3.26) 



We are now in a position to calculate the exponent v* in 1-loop order. The /^-function is 
at this order 




(3(b) = -eb + 6 2 Res + 

= -£& + 26 2 + 0(b 3 ) + C(6 2 ^) , 
and the scaling function u{b*) becomes 

„(*•) = i-iR»(0)£ + 0( 



+ C(6 3 ) + C(6 2 £) 



(3.27) 



1 £ 

2 ~ I' 



Res 



Res 



\ + l + 0(e*) 




+ 0(e 2 ) 



(3.28) 



This renormalization scheme is also used in <y? 4 -theory. At 1-loop order, no renormaliza- 
tion of the wave-function is necessary. Only the reduced "temperature" t is renormalized. 
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There is another scheme, equally useful, to perform the renormalization of polymers, 
which is also used in the broader context of polymerized membranes. This scheme also 
works for infinite membranes. Naturally, for infinite membranes, no renormalization of 
t can occur as it is identically 0. It is also known for the renormalization of standard 
field-theories that one has the choice to work either in a massive (t ^ 0) or a massless 
(t = 0) scheme. 

For the polymer model, let us find a renormalization scheme where t is not renormal- 
ized, and therefore the limit t — > can be taken without problem. The key observation is 
that only the combinations ZZ t and ZbZ d ^ 2 enter the renormalization group calculations, 
and these combinations are left invariant by changing the Z-factors to 

Z' t = l, 

Z' = ZZ t , (3.29) 

r/l ry ry£ — 2 

L h — L h L t 



For a derivation of this property as a consequence of the rescaling-invariance of the un- 
derlying Hamiltonian, see Appendix B. In terms of the modified Z-factors, we obtain 



l + bi-lnZ' + lbi-lnZ' 



v{b) = l --\p>{b)^Z> . (3.30) 

This is the scheme used by David, Duplantier and Guitter [23,24], and by David and Wiese 
[25,26] in the context of polymerized membranes, where it is the most suitable for higher 
loop calculations. On the other hand, it may lead to some confusion as it necessitates 
a renormalization of the field, even in the case of polymers. This may not have been 
expected from the 1 to 1 correspondence on the level of diagrams for the N — > limit 
of <^ 4 -theory, and polymers. As shown above, the two schemes are completely equivalent 
and one may use the one better suited to the problem at hand. 

Let us stress another important difference between the two approaches. This is most 
easily done by using the multilocal operator product expansion (MOPE) introduced in 
Refs. [23,24], and heavily used in Refs. [43,25,26]. To this end, let us write the multilocal 

operator 5 d (r(x) — r(y)) as x * * y . Divergences in the perturbation expansion then 

occur when distances become small. The first such configuration is the contraction of the 

end-points of one dipole, which we shall denote by (^) . A derivation of the contribution 

of this diagram in the more general context of self-avoiding membranes is given in appendix 
C. Specializing to polymers gives the result 



\x-y\ 



d ' 2 l-\\x-y\ 1 - d ' 2 ++... , (3.31) 



2 



where 



+ == ^(Vr) 2 • (3.32) 

The divergence proportional to the operator 1 is subtracted by analytical continuation. 
In the absence of any boundary and for infinite membranes, this term has no effect on 
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the renormalization functions. The second term is more serious and has to be subtracted. 
This is done by renormalization of the field, thus introducing the renormalization factor 



Z = 1 + — . 

2e 



(3.33) 



Upon expanding the Hamiltonian, this yields a counter-term proportional to +, which 
cancels the divergence. 

Let us now study the renormalization of the coupling constant in this scheme. Using 
the MOPE, we can write down the following two UV-divergent configurations 



and 



(3.34) 



from which we shall extract terms proportional to the interaction • • , which we denote 

as 



• • ) , and 



The first is written in the notation of polymer theory as 



■X 



(3.35) 



(3.36) 



Indeed this diagram was subtracted when we renormalized the interaction in Eq. (3.24), 
where we also subtracted the term, 



(3.37) 




The MOPE now tells us that 



x 



(3.38) 



as will also be proved in the context of membranes (see appendix C, or Refs. [23-25], where 
this is discussed in some detail). This result implies that having introduced a counter- 
term for (^T) , i.e. a renormalization of the field, no counter-term for the diagrams in Eq. 
(3.37) is needed. 

We can check for consistency by comparing the /9-functions from the two schemes at 
1-loop order. In the massive scheme, we had 



(3(b) = -eb + b 2 Res ^ + j 




+ 0(b 3 ) + 0{b 2 e) . 



(3.39) 



In the massless scheme, we obtain 
(3(b) = -eb + b 2 (iW !' 



->-2 ReS V. 



+ ) + 0(b 3 ) + 0(b 2 e) . (3.40) 
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It is now easy to see that expressions (3.39) and (3.40) are equivalent up to order 0(b 3 ) 
and 0(b 2 e), since 



Res 



. ) = Res 




and 



--Res 
2 



+ 



Res 



= - A = l + 0{e). 




= 1 



(3.41) 



(3.42) 



ik(r(x)—r(y)) 



Another observation is that in the massless scheme, vertex operators like e 
finite, whereas in the massive scheme they require an additional renormalization. 



are 



4 Generalization to N colors 

Having performed a careful analysis of the different renormalization schemes, we are now 
in a position to generalize to the case N > 0, i.e. to an arbitrary number of self-avoiding 
polymer loops. To this aim, we introduce polymers of iV different colors, and for the 
time-being, work in the massive scheme. In addition to ^ j ^ , which renormalizes the 

chemical potential t, there is now a second contribution, namely 

OO ■ <"> 

This diagram is easily factorized as 

OO = O x O • < 4 2 > 

and is therefore equivalent to the digram already encountered in Eq. (3.12) and absorbed 
in Z t (for N = in Eq. (3.16)). 

Let us now determine the combinatorial factor: A configuration 



OO 



(4.3) 



can be made out of 1 polymer in 2 different ways or out of two polymers. The latter 
comes with an additional factor of N, accounting for the N different colors introduced 
above. Z t is thus modified to 

* = 1 " l Res (O) ( l + T ) + OOT 

= 1 + Y e (l + j) + 0(b 2 ) . (4.4) 

This is indeed the same combinatorial factor as derived from iV-component y9 4 -theory. 
For the renormalization of the coupling-constant, in addition to 




and TV{ , (4-5) 
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there is the possibility that an additional loop mediates the interaction between two given 
polymers, described by a configuration 



The configurations in Eq. (4.5) are realized in 4 different ways each in the high-temperature 
expansion, while for Eq. (4.6) there is only one realization which comes with a factor of 
N for the N different colors. Z& is therefore modified to 




i K8 + A0 

= 1 + i — - • 4 - 7 

Ae 

Evaluating the critical exponent v* as before now yields 

* 1 e2 + N 

It is again possible to switch to the massless scheme. At this stage this is not very 
enlightening, as for polymers all diagrams are essentially equivalent. We will therefore 
discuss this scheme in the context of membranes, which is introduced in the next section. 



5 Generalization to membranes 

We shall now introduce a generalization to polymerized tethered membranes. Formally, 
one generalizes the function r(x) to [12-17,20-26,43] 

r : x G 1R D -> r(x) £ JR d , (5.1) 

with membranes obtained from D = 2. While perturbation theory is always singular for 
D — 2, it is possible to perform an analytical continuation in the inner dimension for 
< D < 2. It is now possible to make an e-expansion, where 

2 - D 

e = 2D — v d , with u = — - — , (5.2) 

about any point (D, d) for which e = (see Fig. 5.1). (The final results are then evaluated 
for D = 2 or 1.) 

We shall re-write the Hamiltonian as 

n = I \(yr(x)f + b^Z b jf jf 5 d (r(x) - r(y)) + tZ t tt , (5.3) 

where the normalization of the integration measure has been chosen for convenience such 
that 

1 f 2tt d / 2 
= -/dx, where S D = fm . (5.4) 
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Figure 5.1: The critical curve e(D, d) = 0. The dashed line corresponds to the standard polymer pertur- 
bation theory, critical in d = 4. 



For radial symmetric functions, the integration is then simply 

Jx JO X 

and the volume of the membrane is 

n = Jd D x = s D J^ . 



(5.5) 



(5.6) 



Note that for D — 1, normalizations agree with those used for polymers in the preceding 
sections. 

The factor 1/(2 — D) is introduced in Eq. (5.3) in order to obtain 

C (x — y) := — ^- (r(x) — r{y)) 2 ^j — \x — y\ 2 ~ D + subleading terms . (5.7) 

To calculate the next to leading term in the free correlator of Eq. (5.7), note that for any 
function Cq(x) defined on a closed compact manifold 



/ AC (:r) 

Jx 



. 



(5.8) 



For the correlator to satisfy the above condition, the usual Laplace equation, AC (x) 
S D (x), has to be modified to 



(5.9) 
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where f2 is the volume of the compact manifold. The numerical prefactors come from 
our choice of normalizations in Eq. (5.3). In the infinite-volume limit, the correction term 
disappears, and the usual equation is regained. It is easy to deduce that 

__/_.\\2\ i \2-d vqSd, ,2 



Cq(x) — — ( -(r(x) — r(y)) ) — \x — y\ t^ft\ x ~ v\ + subleading terms . (5.10) 

a \ 2 / D\l 

The coefficient of the correction term clearly agrees for D = 1 with the exact result for 
closed polymers in Eq. (3.4). 

The considerations of section 3 can now be generalized to the case of membranes. The 
free partition function for a single polymer, i.e. the sum over all sizes of a non-interacting 
polymer, Eq. (3.9), is generalized to 



( 0) = cCD) rdn nn _ uod/D m 
1 d J n 



< D \ r (L-i)f/^ . (s.ii) 



D \D 



We have chosen to integrate over a logarithmic scale, — = 7517. To emphasize the 



X 



arbitrariness of this choice, we have included an additional factor of c(D), which is further 
discussed in the next section. This factor is important, as it also appears in the ratio of 
divergences due to self-interactions of one membrane, and those of interactions with 
other membranes. The factor Q in the integrand of the above equation originates from 
the possible choices of a point xq on the membrane, while the factor 

Q-vod/D „ (s^r(xo))) Q , (5.12) 

is the probability that at this point the membrane is attached to a given point in space. As 
usual, we have introduced a chemical potential proportional to the size of the membrane. 

Let us now generalize Eq. (3.12) for the effect of one 5 d -insertion from the expansion 
of the interaction. For the time being, we fix the size of the membrane to Q, and evaluate 

ix Jy 



5 d (r(x)-r(y))) o . (5.13) 



This integral is (see Eq. (3.7)) 

I- I Co(x)-^ , (5.14) 



Sd j x 

and we have to remove all UV-divergent contributions. To do so, we expand Co(x)~ d >' 2 
for small x. Up to UV-convergent terms, this is (using Eq. (5.10)) 

£ / ^ (x-^ + i^aP-x* + ..) . (5.15) 
S D J x V 2 Dn I 

The first term is strongly UV-divergent and has to be subtracted by a finite part prescrip- 
tion, while the second is (up to terms of order e°) equal to 

-n £/D . (5.16) 
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Note that we have cut off the integraf at the upper bound aj max = VL 1 / . This procedure 
may appear rather crude, but the residue of the pole in 1/e is not affected [43]. 
Upon integrating over all scales, the partition function (to first order) reads 



2 W = <D) 



1 



tt Q-"od/D e -tn 



D J tt 

which upon integration over Q results in 



= Z{ 0) 



l + -{Qt) £ / D + 

€ 



(5.17) 



1 + Ye + - 



(5.18) 



Note the difference in factor of 2 between Eqs. (5.16) and (5.18), which is due to a subtlety 
known to result from nested integrations in standard field theories. This factor of 2 can 
also be interpreted as being geometric. The counter-term is only needed in the half-sector 
x < Vt l l D and not in the half-sector x > Q}/ D . (See also the calculations in Refs. [20,25]). 
Introducing now a counter-term for t yields 



Z t = 1 + — . 

2s 



(5.19) 



The bare and renormalized quantities are now related by generalizing Eq. (3.14) to 



r = VZr , 

— Z t t , 
b = /2 £ Z d / 2 Z b b , 



(5.20) 



leading to the renormalization group functions (compare with Eqs. (3.25) and (3.26)) 



-eb 



l + 6|lnZ b + f6|lnZ ' 



(5.21) 
(5.22) 



(The derivation is given in appendix A.) 

The combinations ZZ^ D ^ D and Z b Z d / 2 , which enter the renormalization group cal- 
culations, are left invariant by changing the Z-factors to 



Z>t — Zt/ Z a , 

z' = zz£- d v d 
zt = z h zf/ D - 2 . 



(5.23) 



For a derivation of this property as a consequence of the rescaling-invariance of the un- 
derlying Hamiltonian, see Appendix B. 

In order to eliminate the renormalization of t, we chose 



Z a — Z, 



t i 



(5.24) 
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resulting in 



z' t = i, 
z' = zz 



(2-D)/D 



(5.25) 



Zl — ZhZ, 



e/D-2 



and the renormalization group functions 

m = 

1/(6) = 



-eb 



l + 6|ln^ + ffe|lnZ' ' 
2 2 1 y J db 



With this change of variables, Eq. (5.19) is replaced by Z' t = 1, and 

2-D& 



Z' = 1 + 



2D e 



(5.26) 



(5.27) 



The above result is precisely that obtained by using the multilocal operator product 
expansion technique for infinite membranes (see appendix D and Ref. [25]), where the 
renormalization factor is calculated from 



Z' = l- (2-D)-Res^^ 



with 

The interpretation of this formula is simple, as 

+ 



Res 



+ 



1 
2~D 



(5.28) 



(5.29) 



is just the diverging contribution form the MOPE, of one <5 d -insertion. For N = 0, the 
renormalization of the coupling constant in the massless scheme is analogously (see section 
3, appendix D, and Ref. [25]) 



Z' b = 1 + -Res^ (i 



(5.30) 



with 



Res( :': 



L_ r JM 

r (s) 



2- D 



(5.31) 



Alternatively, in the massive scheme (Z — 1), 



Z b = ( 1 + -Res(^ 



= 1 + - ( Res 



xZ t 2 



+ 1 



(5.32) 
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Let us now study the generalization to iV components in the massive scheme. Taking 
care of the additional factor of c(D) introduced in the definition of the free partition 
function in Eq. (5.11), Eq. (5.19) is modified in the same manner as Eq. (4.4) to 



Z t = l + 



2e 



1 + 



(5.33) 



There are several possibilities to derive the modification to the renormalization factor of 
the interaction. For a direct derivation generalizing Eq. (4.7) we calculate the diagram 
for membranes, as 



(5.34) 



We already have given the derivation of a similar integral in Eq. (5.11). The only difference 
is that now a second factor Q appears to take into account the additional point which 
moves on the membrane. Integration over Q yields 



o 



c(D) 



D \ D 



r £ l D = c{D)-t £/D + O(e ) . 



(5.35) 



This term again appears with a relative combinatorial factor of N/4 compared to the other 
contributions, as discussed in the polymer case. The renormalization factor Zb therefore 
becomes 



1 + - I Res 



+ 1 + 



jD)N \ 



(5.36) 



It is now easy to derive the renormalization group functions 

c(D)N\ 



(3(b) = -eb + b 2 Res 



+ 1 + 



+ 0(b 3 ) + 0(b 2 e) , (5.37) 



and 



1/(6) = 



2-D 



1 + 



2D 



1 + 



c(D)N" 



0{b 2 ) . 



(5.38) 



At the non-trivial (IR-stable) fixed point, this yields the critical exponent to order e 



2-D 



( 



v 



1 + 



1 + c{D)N/2 



\ 



V 



2D 



Res 



• ) + l + c(D)N/4j 



(5.39) 



representing our central result for the generalized 0(iV)-model, also discussed in the 
introduction in Eq. (1.1). 
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6 The arbitrary factor c(D) 



In calculating the free partition function in Eq. (5.11), we introduced an arbitrary factor 
of c(D). In principle, any function of D which satisfies 

c(l) = 1 , (6.1) 

reproduces the correct result for linear objects. The additional freedom (or ambiguity) 
is apparently a reflection of the non-uniqueness of the generalization to manifolds. Even 
after restricting to the class of hyperspheres, there is a remaining ambiguity in the choice 
of the measure for the size of these manifolds. This arbitrariness carries over to our 
generalization of the O(N) model to iV-colored membranes. (Note also that c(D) is 
independent of the introduction of factors like the -^zd ^ n ^Q- (5-3).) 

We shall focus on two natural choices for c(D). The first possibility is to demand that 
the free partition function in Eq. (5.11), 

s?» = ^>r(£-i) «■/"-', (6.2) 

has the simplest possible form in the sense of depending on only one parameter besides 
t, namely e/D. This implies 

c(D) = D , (6.3) 

which is our favorite choice, and equivalent to the measure d.Q/0 over all scales. 

Another reasonable choice is to demand that the ratio of the probability that one 
membrane touches itself, to the probability that two membranes touch, which is 

A/Nc{D), (6.4) 

be independent of D. This leads to 

c(D) = 1 , (6.5) 

and is equivalent to the measure dx/x over scales. We shall study both choices in the 
next section. 



7 Extrapolations 

This section is devoted to extracting the information about the exponent v in physical 
dimensions, from the general result in Eq. (5.39). (We shall use v, rather than i/*, to 
denote the fixed point value.) As will become apparent, various extrapolation schemes 
are possible, and choosing the best one is almost an art; we shall rely heavily on the 
methods developed in Ref. [25] to which the interested reader is referred to for further 
details and discussion. The general idea is of course to expand about some point (D , d ) 
on the critical curve e(D , d ) = 0. The simplest scheme is to extrapolate towards the 
physical theories for D = 1, 2 and d = 2, 3, . . ., using the expansion parameters D — D 
and d — do- However, as shown in Ref. [25], this set of expansion parameters is not 
optimal, and better results are obtained by using D c (d) = ^ and e(D, d) = 2D — ^ip-d. 
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Furthermore, it is advantageous to make expansions for quantities such as vd or v(d + 2) 
rather than v. 

We shall mention three such extrapolation schemes which are based on the remarkable 
fact that three different expressions for v coincide with v = (2 — D)/2 on the critical line 
e = 2D - dv = 0: 

(1) The mean- field result 

9 D 

(7.1) 



2D 

~~d~ 



which in the context of polymers and membranes is known as the Gaussian varia- 
tional approximation [8,44-46]. Eq. (5.39) can then be re-organized as 

vd = 2D + a(D)e , (7.2) 

leading to an expansion about the mean- field result, which can be plotted as a 
function of the expansion point D. 

(2) The Flory expression 



2 + D 
^ " YTd ' 

can similarly be used as the basis for the expansion of the quantity v(d + 2). 



(7.3) 



(3) The large N limit of the theory can be solved exactly, as will be demonstrated in 
the next section. The corresponding exact result 

D 

Vn-><x> = Tn~ -i (7-4) 



d- 



2D ' 
2-D 



is also equal to v for e = 0, and can be used as a basis for expansion of the quantity 
v(d- 2D ' 



2-D 



After selecting one of these schemes, the next step is to re-express Eq. (5.39) in terms 

or example, fo 

1 + c(D )N/2 



of D c (d) = ^ and e = 2D - ^p-d. For example, for v{d + 2), the final result is 



v{d + 2) = 2 + D c (d) 



Res 



+ 1 + c(D )N/A 



Do 
2 



2 + D, 

2Dn 



e(D,d) . 
(7.5) 
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Figure 7.1: Results of the extrapolations for v for polymers and the Ising model in two dimensions. 
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Figure 7.2: Results of the extrapolations for v for the O(N) model in 3 dimensions. 
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Figure 7.3: Results of the extrapolations for v for the O(N) model in 3 dimensions, linearized in N. 

If we are interested in the field theory (D — 1) in d — 3, we have to evaluate the above 
expression for 5 = 1/2. However, we are still free to choose the expansion point along the 
critical curve, which then fixes D . As the expansion point is varied, different values for 
v(d+2) (i.e. 5v for d — 3) are obtained, as plotted in Fig. 7.4. The criterion for selecting a 
value for v from such curves is that of minimal sensitivity to the expansion point D . We 
thus evaluate v at the extrema of the curves. The broadness of the extremum provides a 
measure of the goodness of the result, and the expansion scheme. The robustness of this 

0.7t 1 



0.65- 




0.55 



0.5 0.5 1 1.5 2 

D 

o 

Figure 7.4: Extrapolations for v from the expansion of v (d + 2) to (D = 1, d = 3), for N = (lower curve) 
and N = 1 (upper curve). The straight lines are a guide for the eye and represent the known results of 
0.59 and 0.63 respectively. 
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Figure 7.5: Extrapolations for j^v(N) from the expansion of vd to D = 1, d = 3. This yields the flattest 
plateau encountered in all extrapolations, and there is no difference to calculations in higher loop order 
(see [2]). 

0.8t 1 




0.5- 



0-4 0.5 1 1.5 2 

Figure 7.6: Extrapolations for v from the expansion of vd to D = 1, d = 3, N = 0, . . . , 3, c(D) = D. The 
dashed lines are marks for the eye and represent the exact values. 

choice in the case of N = was explicitly checked in Ref. [25], by going to the second 
order. (For additional discussions of such "plateau phenomena" see Sec. 12.3 of Ref. [25].) 

While we examined several such curves, only a selection is reproduced in Figs. 7.4, 
7.5, and 7.6. We start by checking the method for polymers (N = 0) and the Ising model 



25 



0.8 




0.5 



0.4 , 



0.5 



D, 



o 



Figure 7.7: Extrapolations for v after linearization in N from the expansion of v(d+ 2) to (D = 1, d = 3) 
for N = 0, . . . , 3, with c(D) = D. The dashed lines are marks for the eye and represent the exact values. 



(N — 1) in d — 2, where the exact values are known [y = 3/4 and 1 respectively.) The 
results are given in Table 7.1. Only the extrapolation for vd yields acceptable results; that 
of v{d + 2) is bad for N — 1, and v is even worse in both cases. Furthermore, we observe 
that c(D) = D gives exponents closer to the correct value. Based on this experience, we 
focus on the expansions for vd in d — 3, with c(D) = D, shown in Fig. 7.6 for N — 0, 1, 
2 and 3. The values of v extracted from the maxima are given in Table 7.2, along with 
their best known estimates from Ref . [2] . Our results are clearly better than the standard 

1- loop expansion of 

1 N + 2 

There are, however, systematic differences in the trends. In particular, we make the 
observation that the "exact" exponents in the range < iV < 3 approximately fall on 
a straight line with slope 0.042 ± 0.003; while our results have a perceptible downward 
curvature. It may be that the 1-loop results are somehow most suited to give the exponents 
at small N. A similar suggestion was made in Ref. [47] in the context of the non-linear 
sigma model. Based on these observations, we now pursue an alternative expansion, that 
searches for the slope of z^(iV) at N = 0. Fig. 7.5 shows the result for du(N)/dN\ N=0 
as a function of the expansion point. A very flat and well defined plateau is obtained, 
with a value of 0.042 in excellent agreement with the slope quoted earlier. Of course, to 
get the absolute value of the exponents, we also need to specify v(N = 0). As discussed 
in Ref. [25], the best expansion quantity for this purpose is v[d + 2) leading to u(N = 
0) = 0.59 at the 1-loop order. (This exponent is also obtained if one demands the 1- and 

2- loop results to be equal.) 
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8 The limit N — > oo and other approximations 

As in the O(N) model, it is possible to derive the dominant behavior for large N exactly. 
In the standard (/? 4 -theory, one starting point is the observation that 

Sr(r)) = (S*(r)) 2 , (8.1) 

since in the limit iV — > oo, spin-components of different colors decouple [2,48,49]. This is 
also known as the random phase approximation (RPA). 

Here, we pursue a slightly different approach, based on the diagrammatic expansion. 
Note that for N — > oo, only simply connected configurations survive. (The vertices 
are made out of membranes, the links out of ^-interactions.) For example, the dia- 
gram which is doubly connected, and the diagram ( j ) which in- 
cludes a self-interaction, each have one factor of iV less than the simply connected graph 
(^^■■^^■■^^ • The fading diagrams for the membrane density at the origin are then 
given by 

f := ( //(rW) ) = + . + ooo + 6-o + < 82 > 

The above sum can be converted into a self-consistent equation for / by noting the 
following: Successive diagrams can be obtained from the first (bare) diagram by adding 
to each point of a manifold a structure that is equivalent to / itself. This is equivalent 
to working with a single non-interacting manifold for which the chemical potential to is 
replaced by an effective value of to + bof. Calculation of / for this manifold proceeds 
exactly as in Eq. (5.11), and results in the integral 

The above integral is strongly UV- divergent, and leads to a form 

f = B(t + b f) 2 -^ d ~ 1 + A, (8.4) 

where 

B = v(l- ^-d) . (8.5) 



2D 

The strong UV- divergence, controlled with an explicit UV-cutoff, is absorbed in the con- 
stant A. It is usually dropped in a dimensional regularization scheme, as in Eq. (5.11). 

The radius of gyration R is now related to / as follows: From Eq. (8.3) we note that 
to + bof is the physical chemical potential conjugate to Q, thus leading to a typical volume 
of Q ~ l/(to+bof). Since there are no self-interactions in the effective manifold introduced 
above, its radius can be related to the volume by R ~ Vt UQ l D . Thus, up to a numerical 
factor which is absorbed into the definition of R, we obtain 

2D 

R-— = to + b f . (8.6) 
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Eliminating / in Eq. (8.4) with the help of Eq. (8.6) yields 

2D , , ID j 

R-— = fa + b A) + b B R^~ d . (8.7) 

Identifying the difference in temperature to the critical theory as 

i=t + b A , (8.8) 

the critical theory is approached upon taking t — > and R — > 00. This occurs if and only 
if d is larger than the lower critical dimension 

2D 

d> d t = . (8.9) 

2-D v 1 

If d is in addition smaller than the upper critical dimension, i.e. 

AD 

d<d u = j-^, (8.10) 

the left hand side of Eq. (8.7) vanishes faster than the independent term on the right 
hand side, and we obtain the scaling relation 

Rrvt- 1 "*-^ . (8.11) 

In the large N limit, the exponent v* is therefore given by 

D 



N^oo ^ 2D 



(8.12) 



2-D 

We can verify that the standard result [2] is correctly reproduced for D — 1 as 



^ooP = l) = j^2 ■ (8-13) 
Note that for d > d u , the leading behavior from Eq. (8.7) is 

R ~ f-T# , (8.14) 

implying the free theory result 

2-D 

= —j- ■ (8-15) 

It is interesting to cast the other approximation schemes introduced in the previous 
section in the language of renormalization factors. The Flory-approximation assumes that 
the elastic energy and the contribution due to self-avoidance scale in the same way. This 
enforces for the renormalization factors Z and Z^ (in the massless scheme with Z t = 1) 
the constraint 

Z = Z b . (8.16) 

Knowing such a relation, the critical exponent v can be calculated explicitly [25]: Suppose 
that 

Z h = Z a , and Z t = 1 . (8.17) 
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From the definition of the /3-function in Eq. (5.21), we obtain 



P(b)§- b ln(Z b Z^) = -e-^. (8.18) 

The second term on the r.h.s. can be neglected upon approaching the critical point b*, 
where the /3-function vanishes. Inserting Eq. (8.17) into Eq. (8.18), solving for P(b)-^ In Z : 
and substituting the result into Eq. (5.22) yields 

v* = — - + . (8.19) 

2 d + 2a K ' 

For the Flory-approximation, a = 1 from Eq. (8.16), and the above expression evaluates 
to 

2 + D , 

^ = YTd ■ (8 ' 20) 

The last approximation scheme that we shall discuss is a mean-field limit. As is well 
known from the mean-field approximation to y? 4 -theory, minimizing an expansion of the 
order parameter (e.g. by a saddle point method) results in an exponent a = describing 
the singularity in the heat capacity. Assuming that the singular part of the free energy 
scales as 

/sing ~ ~ ^- dU ' ,D ~ t dV *' D , (8.21) 

leads to a generalized heat capacity exponent for manifolds given by 

a = 2-^. (8.22) 

A discontinuous (but non-diverging) heat capacity then leads to 

2D 

%f = -j , (8.23) 

which coincides with the result obtained by considering large d and N = [44-46]. In 
such a limit, and in a massless scheme, corrections due to self-avoidance are strongly 
suppressed [25], and we have 

Z b = l , Z t = l, Z 1 . (8.24) 

This is equivalent to a = from Eqs. (8.17) and (8.18), and again yields Eq. (8.23). 



9 Low temperature expansions of the Ising model 

In Sec. 2 we demonstrated that the high temperature expansion of the O(N) spin model 
naturally leads to a sum over iV-colored loops (D — 1); motivating the later generalization 
to manifolds (arbitrary D). For the Ising model (N = 1), a related description can be 
obtained from a low temperature expansion. Excitations to the uniform (up or down 
pointing) ground state are in the from of droplets of spins of opposite sign. The energy 
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cost of each droplet is proportional to its boundary, i.e. again weighted by a Boltzmann 
factor of the form 



Thus a low temperature representation of the d- dimensional Ising partition function is 
obtained by summing over all closed surfaces of dimension D = d — 1. For d = 2, the 
high and low temperature series are similar, indicating the self-dual nature of the model. 
For d = 3, the low temperature description is a sum over surfaces, which is also the 
high temperature expansion of an Ising lattice gauge theory [11], establishing the duality 
between these two models. 

The non-trivial question is regarding the type of surfaces which dominate the above 
sum. There is certainly no constraint on the internal metric, in contrast to the tethered 
surfaces in D = 2 which have a flat metric. Since the sum includes droplets of all shapes, 
it may be more appropriate to examine fluid membranes. However, there is currently 
no practical scheme for treating interacting fluid membranes, and the excluded volume 
interactions between the membranes are an essential ingredient to avoid overcounting 
configurations. We shall argue that, at least in low dimensions, the sum over tethered 
membranes captures the appropriate physics of the problem. This may appear quite 
surprising at first glance, as for N = 0, surfaces generated from plaquettes on a lattice 
are very different from tethered surfaces. The former are dominated by configurations 
that resemble branched polymers. The large entropy gain of branches is responsible for 
this, and appears as an instability towards formation of spikes in a string theory [50]. 
However, it is possible that for iV > 0, the above instability is replaced by a string of 
bubbles. (Reminiscent of the Raleigh instability of a stream in hydrodynamics.) The 
collection of bubbles is then satisfactorily described by a set of fluctuating hyperspherical 
(tethered) manifolds which is the basic ingredient of our model. The appropriate question 
may be whether tethered membranes sweep out phase space, i.e. form a complete set of 
basis functions in the configuration space of our problem. 

Another issue is whether the sum may be restricted to spheres, or if objects of other 
topologies must also be included. We argue in Appendix E that the dominant contribution 
(and the only one included in perturbation theory) is the one from spheres. 

We shall now test the validity of the above conjecture. As discussed before, our 
generalization to the sum over manifolds is defined only up to the factor c(D). For the 
remainder of this section we make the "least number of variables" choice c(D) = D, for 
the partition function. Singularities of the partition function are characterized by the 
critical exponent a(D,d, N), through 



The equality of the singularities on approaching the critical point from the low or high 
temperature sides, requires 



t-t c 



2-a 




'singular 



a(l, d, 1) = a(d — 1, d, 1) . 



(9.2) 



From the scaling-relation 



a(D,d,l) = 2 



v(D,d,l)d 
D 



(9.3) 
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Figure 9.1: Test of the relation in Eq. (9.4) for the Ising-model in d = 2.25, d = 2.5, d = 2.75 and d = 3. 
The upper curves are from the high temperature representations of the Ising model (D = 1), while the 
lower curves are from the low temperature expansion (D — d— 1) as explained in the text. The exponent 
v is extrapolated from the maximum of each curve. These curves are based on the scheme in which we 
extrapolate vd, linearize in N, and to divide the result by dD. This choice is determined through the 
quality checks of section 7. 

we obtain 

u(l,d,l) = ^ju(d-l,d,l) . (9.4) 

A more general identity is obtained by considering a generalized class of Ising models, 
whose interactions are defined on .D-dimensional primitive elements of a lattice [11]. In 
the standard Ising model the interactions occur on bonds (D = 1), while in Ising lattice 
gauge theory interactions are placed on plaquettes (D = 2). Equating the singularities 
from the high and low temperature expansions of such models, and assuming a single 
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Figure 9.2: Extrapolations for the Ising model and its dual in D 
extrapolation variables 1 /D and d. 



1, d = 3, using c(D) = D, and as 



continuous phase transition, yields 

u(D,d,l) = j^u(d- D,d,l) . (9.5) 

The conjectured identity in Eq. (9.4) was tested numerically, and the results are pre- 
sented in Fig. 9.1. The extrapolated exponents (the maxima of the curves) from the dual 
high and low temperature expansions are in excellent agreement. Indeed, one could hardly 
expect better from a 1-loop calculation. Nevertheless, higher-loop calculations would be 
useful to check this surprising hypothesis. One of the peculiarities of the 1-loop extrapo- 
lations presented in Fig. 9.1 is that any crossing of the curves from the dual models occurs 
at the mean-field value of 2D/d from Eq. (8.23). This is accidental, but present for all 
of our "good" extrapolation variables. An explicit counter example is given on Fig. 9.2, 
where we used the "bad" extrapolation- variables 1/D and d. They are bad, as there is 
no pronounced maximum to estimate v. Nevertheless, the intersection of the two dual 
models occurs at v = 0.6, not too far from the exact result of 0.6315. 



10 Cubic anisotropy 

Up to now, we assumed that the interaction between manifolds is independent of their 
color. This is a consequence of the rotational invariance of the underlying spin model 
introduced in Sec. 2. This equality of interactions does not have to hold in a system of 
polymer loops. In the context of the (^ 4 -theory, unequal interactions result from the break- 
ing of rotational symmetry, e.g. through the introduction of cubic anisotropy, as discussed 
in Refs. [49,51-53]. In the microscopic spin model of Sec. 2, the independence of the in- 
teractions between loops from their colors emerges as a consequence of the normalization 
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Figure 10.1: Regions of stability for the bare model with cubic anisotropy (the white portions are stable): 
The bare model is unstable if u + Nb < for u > 0, and if u + b < for u < 0. Note that this diagram 
changes upon renormalization. 



condition 



E^ 2 = i. 



(10.1) 



If we replace this by the constraint 



(10.2) 



with a ^ 2, which breaks rotational symmetry, this is no longer the case. 

The model in the absence of full spherical symmetry is described by two interaction 
parameters. In addition to b which indicates the interaction between any two membranes 
irrespective of their color, there is a new anisotropic coupling constant u which acts only 
between membranes of the same color. Physically, not all combinations of u and b are 
allowed, as some of them induce a collapse of the membranes. The following two cases 
can be distinguished: 

(i) If u is negative, for a single membrane to avoid itself and not to collapse, the 
condition 



has to hold. This implies that b is positive, and the repulsive interactions between any 
pair of membranes ensures mutual stability. The same condition is obtained in the O(N)- 
model by requiring the stability of the minimum energy state [49]. 

(ii) If u is positive, the stability condition in the 0(iV)-model is [49] 



This places a lower bound on b which is, for more than one color (JV > 1), more restrictive 
than Eq. (10.3), but still admits negative values for b (see Fig. 10.1). For b < 0, mem- 
branes of different color attract; in the extreme limit becoming glued together to form 



u + b > 



(10.3) 



u + Nb > . 



(10.4) 
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a "super- membrane" out of N differently colored membranes. In this limit, the theory 
reduces to an Ising-like system, where the effective number of colors is one. (The corre- 
sponding RG-flows, as discussed below, do indeed tend to an Ising fixed point.) For this 
"super- membrane" not to collapse, we again find the condition in Eq. (10.4). However, 
this picture is quite schematic, and real physical systems are governed by many more 
parameters, and may well behave differently. Studies in polymers [54,55] indicate that 
the precise competition of the attractive and repulsive parts of the interaction potential 
plays a crucial role, sometimes leading to non-universal behavior. 

The above stability arguments are based on energetic considerations, and are expected 
to be modified upon the inclusion of fluctuations, say through a renormalization group 
procedure. In the studies of critical phenomena, a well-known example is the Coleman- 
Weinberg mechanism [49,56], where the RG flows take an apparently stable combination 
of b and u into an unstable regime, indicating that fluctuations destabilize the system. 
In the flow diagrams described below, we also find the reverse behavior in which an 
apparently unstable combination of b and u flows to a stable fixed point. We interpret this 
behavior as indicating that fluctuations actually stabilize the model, a reverse Coleman- 
Weinberg effect, which to our knowledge has not been discussed before. To decide whether 
a system with a given combination of couplings is stable, we first follow the renormalization 
group flow. As perturbative renormalization does not say anything about the strong 
coupling regime, we regard run-away-trajectories in the flow as indicating an unphysical 
situation. If, on the other hand the renormalization group flow tends to a finite and 
completely IR-stable fixed point, we use the "classic" stability analysis discussed above, 
since renormalization has eliminated all fluctuations. Using this criterion we have shaded 
in grey the unphysical regions in the following flow diagrams. 

The derivation of the generalized renormalization group functions is most easily carried 
out in a mixed scheme, in which we absorb contributions due to self-avoidance into Z, 
and those due to interactions with another membrane (proportional to N) into Z t . For 
the original model, this scheme leads to the renormalization factors 

w 2e 2 

2 - Db , 

Z{h) = l + ^W~~e- (1 °- 5) 

The presence of the additional interaction u between membranes of the same color, mod- 
ifies the above result to 

z(M) = i + H^±!±. (io.6) 

To derive the renormalization of the coupling-constants, we note that with this choice 
of Z and Z t , we have eliminated all divergent configurations which in the polymer picture 
and in the MOPE are denoted by 
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respectively. The diagrams that renormalize the interactions are of two classes 



and 



XX' 



(10.8) 



In calculating the contribution to the renormalization of b, the inner loop of the latter 
diagram contributes a factor of N if both interactions are b, and 1 for the two cases when 
one of the interactions is u, resulting in 



b = {i~ £ Z~i fl - ^Res 



b : u — - — I 6 . 



As 



2e 



10.9) 



(Note that because of the choice of the mixed RG scheme, the above relation does not 
reduce for u = to the corresponding one derived earlier.) For the anisotropic coupling 
u, the first class of diagrams can be constructed from two interactions u (1 way), or one 
u and one b (two ways). The second class gives a single contribution proportional to u 2 , 
for the overall result 



u = fi £ Z 2x1 Res 



c(D)\ 
u— u . 



(10.10) 



The renormalized parameters yield the /3-functions at 1-loop order, as 



(3 b (b,u) = -eb + b 2 Res 



+ 1 + 



(10.11) 



and 



P u (b : u) = —eu + u 2 ^Res 
+ bu ( 1 + 2Res( 



+ 1 + 



c(DY 



(10.12) 



Finally, the exponent v is obtained from Eq. (10.6) as 
u(b, u 



2 \ 2D 



10.13) 



For D — > 1, these equations reduce to the renormalization-group functions reproduced 
in Amit [49]. (In [49] there is an additional factor of 2/3 due to the choice of numerical 
constants.) 

As in their standard counterpart, these flow-equations admit 4 fixed points: 
(1) The Gaussian fixed point 

6^ = 0, ^ = 0, (10.14) 

which is always unstable below the upper critical line d u (D). (The following discus- 
sions of stability all pertain to this region.) 
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(2) The Heisenberg fixed point 



b * = 1 

" 1 + M^ + Res((< 



, u* H = 0, (10.15) 



is always stable along the u = axis, and is completely stable as long as 

^e<Re S (@3^|„) ei . (10.16) 

This certainly applies to single polymers and membranes with N = 0, but is also 
the case for the O(N) model as long as iV < 4. 

(3) The Ising fixed point 

b} = , u} = 

is stable if 



;io.i7) 



— )<^- do-*) 

For the standard 0(iV)-model with D = 1 the Ising fixed point is unstable. However, 
since for large d C) the left hand side of the above inequality decays rapidly as 2~ d °/ 2 , 
the Ising fixed point is stable if the expansion point is sufficiently close to D — 2. 

(4) The cubic fixed point is located at 

Kg) - g) 

hc ~ (l + ^ + A(£>)) - ^) + (l + ^) (^P - A(£>)) ' ' 

. = (^-A(D)) 

Uc (l + ^ + A(D)) (A(D) - ^1) + (l + ^) _ ' ' 

where A(D) = Res( 



(10.19) 

The stability of this point under RG flows depends on the parameters D and N, and 
this dependence is not simple. The physical stability of this fixed point according 
to the criteria of Eqs. (10.3) and (10.4) should also be verified. Interestingly, we 
find (at least at 1-loop order) that whenever the fixed point is IR-stable, it falls 
in the physically stable region. Note that there are combinations of D and N, for 
which the cubic fixed point is at infinity in the 1-loop approximation. An explicit 
example is for D — 1 and iV = 0, a scenario relevant to the random bond Ising 
model discussed in the next section. 

The two RG equations admit six different flow patterns, the last four of which do not 
occur in the standard field theory. The domains of the (N, D) plane corresponding to 
each of the following scenarios is plotted in Fig. 10.2: 
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Figure 10.2: Regions with different RG flow patterns in the (N, Z))-plane, as discussed in the text. The 
domain-boundaries are obtained from c(D) = AA{D) for the line separating (iii) and (iv) as well as (ii) 
and (vi); Nc(D) = AA(D) for the line separating (i) and (ii) as well as (iv) and (v). A third boundary is 
given by the vanishing of the denominator in Eq. (10.19). 

(i) For D — 1, i.e. in the case of the O(N) model, and for < N < 4, only the 
Heisenberg fixed point is stable, as indicated in Fig. 10.3. This is the fixed point 
that is usually studied in the context of critical phenomena. 

(ii) For D = 1 and iV > 4, the Heisenberg fixed point is unstable and the system is 
governed by the cubic fixed point (Fig. 10.4). 

(iii) An interesting phase diagram is obtained for N = and 1 < D < 1.29. Then, 
as shown in Fig. 10.5, the Heisenberg and cubic fixed points are both stable, their 
domains of attraction being separated by the axis 6 = 0. 

(iv) Another phase diagram with two completely stable fixed points is obtained for N = 
and D close to 2. Then, as shown in Fig. 10.6, the Heisenberg and Ising fixed points 
are both stable, and there is a phase separatrix passing through the Gaussian and 
the cubic fixed points. 

(v) For D close to 2 and N > N c , where N c vanishes exponentially for D — > 2, the cubic 
fixed point is in the lower right sector. Both the cubic and Ising fixed points are 
stable as indicated in Fig. 10.7. 

(vi) Yet another possibility is that both the Heisenberg and the cubic fixed points are 
unstable, as in Fig. 10.8. This is the case for D close to 2 and N large. Then only 
the Ising fixed point is attractive and controls the critical behavior. 

We may inquire as to how the above flow diagrams, with different stable points and 
stability regions, can occur by continuously moving around in the (N, D) plane. To 
demonstrate this, let us examine the sequence of flow diagrams for iV = and differing 
D. The appropriate flow diagram for D < 1 is that of Fig. 10.3, with the cubic fixed 
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Figure 10.3: RG-flow from Eqs. (10.11) and (10.12) in domain (i), e.g. D = 1 and N < 4; shaded regions 
are unstable. 




Figure 10.4: RG-flow from Eqs. (10.11) and (10.12) in domain (ii), e.g. D = 1 and N > 4; shaded regions 
are unstable. 



point going to infinity as D — > 1. For D > 1 this fixed point reappears in the upper left 
sector, as in Fig. 10.5. Along the way, it coincides with the stable fixed point at infinity, 
and they exchange stability. This mechanism of changing stability from one fixed point 
to another is quite general, and occurs again when the cubic and Ising fixed points merge 
at D pa 1.29. For 1.29 < D < 2, the appropriate flow diagram is that of Fig. 10.6. The 
cubic fixed point continues to approach the Heisenberg one as D — > 2, resulting in very 
slow approaches to the fixed points in this limit. 
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Figure 10.5: RG-flow from Eqs. (10.11) and (10.12) in domain (iii), e.g. N = and 1 < D < 1.29, shaded 
regions being unstable. 




Figure 10.6: RG flow from Eqs. (10.11) and (10.12) in domain (iv), e.g. N — and D close to 2; shaded 
regions are unstable. 



We may ask whether the expanded picture presented here provides any new insight 
into the behavior of tethered self-avoiding membranes (D — 2, N — 0). The perturbative 
expansion predicts a crumpled phase with an exponent of v — 0.85 < 1 in d — 3 [25,26]. 
Yet many simulations of this system using models of beads and springs [27-29,36,57] seem 
to suggest a flat phase with v — 1. It is thus important to ask if there are additions to the 
standard description of such membranes in Eq. (5.3), which can lead to a flat configuration. 
The most natural candidate is a bending rigidity k, which is automatically generated in 
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Figure 10.7: RG-flow from Eqs. (10.11) and (10.12) in domain (v), e.g. N = 0.2 and D = 1.8, shaded 
regions being unstable. 




Figure 10.8: RG flow from Eqs. (10.11) and (10.12) in domain (vi), e.g. large N and D close to 2; shaded 
regions are unstable. 

models of strings and beads as pointed out in Ref. [58]. However, it is expected that a 
finite k c is needed in a flat phase, while the absence of a crumpled phase in simulations 
may suggest k c = 0. Our modified Hamiltonian with cubic anisotropy indicates that the 
presence of even small u for D — > 2 places the system within the domain of attraction of 
the Ising fixed point in Fig. 10.6. (The separatrix through the cubic fixed point approaches 
the horizontal line as D — > 2.) The exponents that we calculate for the Ising fixed point 
(v > 1 in d — 3, 4 and v = 0.9 in d = 5, v = 0.8 in d — 6, and v — 0.6 in d — 8) are 
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tantalizingly close to those found in the simulations of Grest [33]. Yet it is hard to justify 
the inclusion of a finite u, which is meaningless for a single membrane at N — 0. While 
the presence of additional membranes does limit the bending of the membranes around 
it, the net effect is much more than just a simple bending rigidity, as related constraints 
appear on all length scales. 



11 The random bond Ising model 

In this section we analyze in greater detail the model for N = 0. The N — > limit is 
interesting, not only because of its relevance to self-avoiding polymers and membranes, but 
also for its relation to the Ising model with bond disorder. To show the latter connection, 
we start with the field theory description of the random bond Ising model, with the 
Hamiltonian 

"I ,2 1 



H 



2 (VS(r)y + -(t + rj(r))S'(r)+uS%r) 



(11.1) 



where r](r) is a quenched random variable (with rj(r) = 0). Expectation values with 
quenched disorder can be calculated from a partition function that is replicated N times, 
in the limit N — > (for a review see e.g. [59]). Averaging the replicated weight over 
the Gaussian random variable r](r), with r](r)r](r') = 2cr5 d (r — r'), induces an interaction 
between different replicas with an effective Hamiltonian 

HN = I £ \\ (yS a {r)f + tl%(r) + <(r)| - a £ ^(r)^(r) . (11.2) 

The replicated system is thus controlled by a Hamiltonian with positive cubic anisotropy 
u, but negative b = —a. 

A key result in the study of random bond systems is the "Harris criterion" [60], which 
states that randomness is relevant as long as the heat capacity exponent a is positive. 
This is the case for the Ising model, and therefore new critical behavior is expected for 
the random bond system. In the usual field theory treatments [59,61-63], there is no 
fixed point at the 1-loop order. This is due to the vanishing of the denominator in Eq. 
(10.19) for D — 1. However, we now have the option of searching for a stable fixed point 
by expanding about D ^ 1. Indeed, for = and 1 < D < 1.29, the cubic fixed point 
lies in the upper left sector (u > and b < 0) and is completely stable, see Figs. 10.5 and 
11.1. The extrapolation for v at the cubic fixed point is plotted in Fig. 11.2, where it is 
compared to the results for the Heisenberg and Ising fixed points. The divergence of is, 
upon approaching D = 1 from above, is a result of the cubic fixed point going to infinity 
as mentioned earlier. Upon increasing D, the Ising and cubic fixed points approach each 
other, and merge for D = 1.29. For larger values of D, the cubic fixed point is to the 
right of the Ising fixed point (b* > 0) and only the latter is stable. Given this structure, 
no plateau can be found for a numerical estimate of the random bond exponent z/doj an d 
we can only posit the inequality 

^DO > Rising • (11-3) 

While Eq. (11.3) was derived at 1-loop order, it should hold at higher orders, if no drastic 
subleading corrections appear; since it merely depends on the general structure of the 
renormalization group flow. 
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Figure 11.1: Renormalization group flows for D = 1.1 and N = 0. The cubic fixed point is marked by a 
diamond 

Higher loop calculations of the random bond Ising model can be used to expand the 
exponent v in powers of \fe. Eq. (11.3) can then be compared to the 2- loop result [63], 
which in d — 3 reads 



This is slightly larger than the corresponding Ising value at 1-loop order (y ~ 0.583), 
but smaller than the best known Ising result [y ~ 0.631). Three- and four-loop calcula- 
tions were performed in Refs. [64] and [65] respectively; the latter gives v = 0.6714 and 
v = 0.6680 depending on the resummation-method used. Two-loop calculations in fixed 
dimension d = 3 [66] yield v = 0.678. A similar result is obtained within a modified 
Pade-Borel approximation of 3-loop results in Ref. [67], yielding v = 0.666. We also note 
that Ref. [67] stresses the existence of very slow transients in the renormalization group 
flow for small e. The extremely retarded crossover to the random bond fixed point may 
quantitatively explain the Monte-Carlo data which seem to suggest a disorder-dependent, 
and thus non-universal, behavior for the critical exponents. The same characteristic flow 
is also present at 1-loop order for our random bond fixed point candidate, e.g. for D — 1.1 
as in Fig. 11.1, further validating our approach. Within error-bars, all these results are 
consistent with a = 0, i.e. the border-line value of the Harris criterion [60], corresponding 





to 



2 



(11.5) 
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Figure 11.2: Extrapolations of v from the expansion of vd with c(D) = D, for the 0(7V)-model in three 
dimensions. The values of the Heisenberg fixed point for N = are compared to those of the Ising and 
cubic fixed points. The crossing of the latter two curves yields an estimate of v = 0.6315 for the d = 3 
Ising model. 



There is also an exact result [68] that the exponent v in any random system must be 
greater than or equal to 2/d. Generalizing the "Harris criterion" bound to manifolds 
yields a limiting value of 2D/d, which is also the mean-field value discussed earlier. 

The duality between the high and low temperature expansions of the Ising model 
remains valid in the presence of random bonds. As discussed above, the high temperature 
expansion can be presented geometrically as a theory of loops with self-avoidance (u > 0), 
and mutual attraction (b < 0). We can develop a related low temperature expansion as 
follows: Starting with an ordered ground state, the partition function Z[{J}] can be 
expanded as a sum over contributions of droplets of the opposite spin; each element of 
the droplet surface crossing a local random bond J makes a contribution of exp(— 2(3 J). 
A replicated description is obtained for ^[{7}]^ as a sum over droplets of iV different 
colors. The next step is to average over all the random bonds {J}: A specific bond may 
be crossed m = 0, 1, . . . , N times for a given term in Z[{ J}] N . Assuming that the random 
bonds are independently chosen from a Gaussian distribution of width a, each random 
bond contributes a factor of 



exp(— 2(3mJ) = exp 



-2(3 ( J -(3<j)m + 4(3 2 a 



(11.6) 



The first term in the exponent on the right hand side can be regarded as a shift in the bond 
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Figure 11.3: Extrapolations of v from the expansion of vd with c(D) = D, for the random bond Ising- 
modcl in d = 3, from the high and low temperature descriptions. Note that the intersection of the curves 
(from these dual descriptions) at the mean-field value of 2/3 is a property of the 1-loop expansion, as 
discussed in section 9. 



energy due to randomness. The second term represents a pairwise attraction between the 
m manifolds of different colors. Thus the quench averaged low temperature description is 
of a set of D = d — 1 dimensional self-avoiding droplets, with mutual attractions between 
droplets of different colors. (This is easily generalized to models with interactions defined 
on other elementary manifolds.) 

We can now make the conjecture that the low temperature sums are not drastically 
modified by restricting the droplets to tethered membranes. This will again lead to 
the exponent identity in Eq. (9.5). The extrapolations from the dual descriptions of the 
random bond Ising model (at the cubic fixed point) are presented in Fig. 11.3. Because of 
the absence of a plateau, there are no clear points where these curves can be compared. 
The intersection of the two curves provides a specific point for extracting exponents of the 
random bond model. However, as discussed in section 9, at the 1-loop level this occurs 
at the mean-field value of 2/3. 

In the previous sections we relied on a plateau in the extrapolation curves to obtain 
numerical values of the exponents. Another method for selecting a specific point of these 
curves is to look at intersection points. For example, in Fig. 11.2 the curves corresponding 
to the cubic and Ising fixed points intersect at D — 1.29. Since the intersection occurs 
when the two fixed points coalesce on the 6 = axis, they both have Ising symmetry at 
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this junction. We may hope that this point yields precise exponents as all ambiguity in 
the expansion point is removed. The actual numerical prediction of 



v* = 0.6315, (11.7) 

is indistinguishable from higher-order calculations [2]. However, the method is not insen- 
sitive to the extrapolation variables used, and it is thus unclear if more precise results are 
obtained in higher order calculations. 

Finally, we note that the above considerations are easily generalized to multi-component 
spins subject to random bonds. We shall denote the number of components of the field 
by p, reserving the symbol N for the number of replicas, as in the random bond Ising 
model (corresponding to p — 1). Starting as in Eq. (11.1) from 



K p = J E 2 ( V ^W) 2 + 2 (t + ^ (r)) ^ (r)2 + U £ WW , (H.8) 

r i=l L J i,j=l 

yields in analogy with Eq. (11.2), after replicating and averaging over r], 

r N v M p f -i N p N v 

Jr a=li=l LZ z J a=li,j=l a,/3=li,j=l 



' j( rf 



(H-9) 

In the resulting O(N) x 0(p)-model, rotational symmetry is broken in the 0(iV)-sector, 
but remains intact in the 0(j9)-sector. We can next develop a geometrical description 
of the high temperature expansion of this field theory in the language of self-avoiding 
polymer loops, which are then generalized to membranes. It is then easy to see that in 
a perturbative expansion, each closed loop contributes an additional factor of p (as in 
the standard 0(iV)-model, where every closed loop is accompanied by a factor of N). 
All previous results are thus simply generalized by replacing c(D) in the renormalization 
group expressions with p x c(D). Since, according to the Harris criteorion [60], non-trivial 
random bond exponents are obtained in the physical dimensions of d = 2 and d — 3, only 
for p < 2, we shall not pursue this analogy further. 
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A Derivation of the RG-equations 

In this section, we give a derivation of the renormalization group functions. Starting from 

b = b Z^Z- d / 2 fi- £ , (A.l) 
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the /9-function is given through the variation of the renormalized coupling, at fixed values 
of the bare coupling and bare chemical potential, as 



(3(b) := 11 



d_ 
d/i 



b . 



(A.2) 



From the derivative of Eq. (A.l) with respect to /i, we obtain 

(3(b) (l + b^- b ln(Z b Z d / 2 )^ = -eb , 



(A.3) 



which solving for (3(b), yields 



-eb 



l + bpn(Z b Z d /*) ■ 



(A.4) 



The scaling exponent v relates the chemical potential t and radius of gyration R 
through 

R ~ r u/D . (A.5) 

To obtain v, we first observe that the dimensionless combination R 2 t^~ , is a function of 
b and t/fi D only, i.e. 

RH^=f(b,t/n D ). (A.6) 
Since in addition, t can be expressed as a function of to and b only, Eq. (A.6) implies that 



rH— 



r, 2-D 



(A.7) 



Next observe that Rfoo D is independent of the renormalization scale fx. Replacing bare 
by renormalized quantities, we obtain a relation for the total derivative with respect to 
/i, as 



= . 



Combining the latter relation with Eq. (A.7) gives 



r, 2-D 



. 



(A.8) 



(A.9) 



The scaling function of the field, describing the scaling of the membrane at the critical 
point with (3(b) — > 0, is thus 



(A.10) 



(Note that the last term cannot be dropped, as the vanishing of the /5-function is canceled 
by the divergence of the derivatives of the Z-factors at this point.) 
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B Reparametrization invariance 



In this appendix, we shall explore the consequences of a reparametrization 

x — > x' = x Z~ 1/D , (B.l) 

on the Hamiltonian 

n = j x \{Vr{x)) 2 + b^Z b jf J S d (r(x) - r(y)) + tZ t Q , (B.2) 

for a self-avoiding membrane (D — 1 for polymers). (The notation for the rescaling 
factor anticipates renormalization factors that we shall introduce next.) The Hamiltonian 
in Eq. (B.2) is in fact not invariant under this rescaling because of the cutoff implicit 
in the interaction. In order to achieve scale invariance, the cutoff, or equivalently the 
renormalization scale /x, must also be rescaled to 

^H' = VZ]J D . (B.3) 

The Hamiltonian then changes to 

2-D 

77 D r 1 r r 

H = yf^ I - 2 (yr\x)Y + b^Z b Z^J D - 2 I J y ~5 d (r(x) - r{y)) + tZ t Z~^ . (B.4) 

Comparing to the original Hamiltonian then identifies the new renormalization group 
factors 

2-D 
Z' = ZZ a ° 

Z[ = ZtZ- 1 (B.5) 

— ^b^t 

As discussed in the main text, see Eqs. (3.29) and (5.25), the renormalization group- 
functions are left unchanged by the transformations in Eqs. (B.5); the most useful case 
being Z a = Z t . 

C Structure of the divergences and the MOPE 

In this appendix we present a more intuitive description of the structure of the divergences, 
and the multilocal operator product expansion (MOPE) used to prove renormalizability 
[23,24] and for explicit calculations [26,25]. This presentation already appears in French 
in Ref . [69] , but is reproduced here for completeness, and convenience of readers. 
We first remark that with our choice of normalizations, the free propagator 

^(^(r(xi)-r(x 2 )) 2 ) o = |x 1 -x 2 | 2 - D , (C.l) 
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is the Coulomb potential in D dimensions. This analogy with electrostatics will help us 
to analyze the structure of the divergences. The interaction part of the Hamiltonian H is 
reminiscent of a dipole, and can be written as 



Wint = 6/ / ~S d (r(x 1 )-r(x 2 )) 

Jx\ J X2 

= b [ [ [ e ik (r(xi)-r(x2)) _ (C.2) 

The next step is to analyze the divergences appearing in the perturbative calculation of 
expectation values of observables. To simplify the calculations, we focus on the partition 
function 

2= E e-« = (e-«-t) . (C.3) 



' o 

all states 



To exhibit the similarity to Coulomb systems, consider the second order term 
-(Hl sk ) = —J J J J J J (e ik{r{xi) - r(x2)) e ip{r{yi) ~ r(y2)) ) 

z J xi J X2 Jyi Jyi J k Jp 

where E c is the Coulomb-energy of a configuration of dipoles with charges ±k, and ±p, 
respectively. More generally, for any Gaussian measure we have 

iEiM«<)\ =e -3Ei lj *i*i(W*i)-^)) a > l)i ( C5 ) 



e 

/ 

Since for any configuration of dipoles, specified by their coordinates and charges, the total 
charge is zero, the Coulomb-energy is bounded from below, i.e. 

E c > . (C.6) 

This implies that 

e~ Ec < 1 , (C.7) 

and that the configurations which contribute most are those with minimal charge. 

Let us apply the above observation to evaluating the integrals in Eq. (C.4). The 
basic idea is to look for classes of configurations which are similar. The integral over the 
parameter which indexes such configurations is the product of a divergent factor, and a 
"representative" operator. For the case of two dipoles, one with charge k and the other 
with charge p — k contracted together, one only sees a simple dipole with charge p from 
far away, i.e. 



, p x e -k 2 (\s\ 2 - D + \t\ 2 - D ) 



(at 



The second factor on the r.h.s. contains the dominant part of the Coulomb energy E c = 
^ 2 (l s | 2 D + \A 2 ~ D ) °f tli e interaction between the two dipoles; s and t are the distances 
between the contracted ends. The integral over k is now factorized, and we obtain 



J k 
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We define the MOPE coefficient, as 



= (| s |2-D + | t |2-^ 



-d/2 



(CIO) 



The MOPE therefore gives a convenient and powerful tool to calculate the dominant and 
all subdominant contributions from singular configurations. 

For the sake of completeness, let us still calculate the two other MOPE coefficients 



used in the text. The first is for small \x — y\, (u 



2-D 



— I ■ e ikr(x) . . e ~ikr(y) . 
Jk 

= f . e ik(r(x)-r(y)) . Q -k 2 \x-y\ 2 -o 
Jk 



= \x-y\ 



l~2 



2d 



-l/gd 



■ [r(x) - r(y)] 2 : + 



3 -k*\x-y\ 2v o 



1 - - : 



[x — y)Vr 



\x-y 



-vod 



1 —\x — y 

2D 1 y 



D—VQd 



x + y 



+ +.... 



(C.11) 



(The normal ordered operator : O : indicates that all self-contractions of O have been 
subtracted.) We also have to show the factorization property of . (^) , which is the 
product of two (^-interactions 



. e ikr(u) . . g-ifcr(x) . . Q ipr(y) . . g - 



z) : . (C.12) 

u xy z 

We want to study the contraction of x, y, and z, and look for all contributions which are 
proportional to 

. . _ J . e ikr(u) . . e -ikr((x+y+z)/3) . ^ (C.13) 

Jk 

The key-observation is that in Eq. (C.12) no contraction with : e ~ tkr ( x ) : contributes. All 
such contractions yield a factor of k, which after integration over k results in derivatives of 
the (^-distribution. This is equivalent to stating that as long as contributions proportional 
to • • are studied, the following factorization property holds: 



(C.14) 



arc 



This is the reason why, in the massless scheme, divergences proportional to (^z 
already eliminated through a counter-term for ( ), i.e. by the renormalization factor Z. 
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D Renormalization for infinite membranes 



In this appendix, we give a short summary of the renormalization procedure for infinite 
membranes, and derive the one-loop counter-terms. This is a simplified version of the 
corresponding section in Ref. [25], which is included here again mainly for completeness. 

Let us start with a single dipole: When its end-points (x, y) are contracted towards a 
point (taken here to be the center-of-mass z = (x + y)/2), the MOPE is 



11 + 



x v---/ y i x 

The first MOPE-coefficients are given explicitly by 



+ + + 



(D.l) 



x V./ y 



1 

+ 

and where -f- denotes the local operator 



\x-y 



—vod 



\x — y 

2D 1 y 



D—i/qcL 



+ = -(Vr(x)) 2 . 



(D.2) 



(D.3) 



-f- ] is, at e — 0, logarithmically 



The integral over the relative distance x — y for 
divergent. 

The simplest contraction resulting in a dipole is when two dipoles coalesce. The 
corresponding MOPE-coefncient is 



(\x\ 2 »° + \y\ 2v ») 



-d/2 



(D.4) 



where x and y are now the relative distances inside the two subsets. Another possibility 

is to consider the contraction . (^) . But as we have shown in appendix C, this does 

not give a new divergence. 

In the next step, counter-terms are introduced to subtract these divergences. We 
have to distinguish between counter-terms for relevant operators and those for marginal 
operators. The former can be defined by analytic continuation, while the latter require a 
subtraction scale. Indeed, the divergence for 1 is given by the integral 



A- 1 <|x-j/|<L 



L dx 



-x 



D—vod 



A- 1 X 



(D.5) 



- ih ( AD - - L '- D ) ■ 



where A is a high-momentum UV-regulator, and L a large distance regulator. For e ~ 0, 
this is UV-divergent but IR-convergent. The simplest way to subtract this divergence is 
therefore to replace the dipole operator by 



(D.6) 
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where • • = la; 



y\ u ° d . This amounts to adding to the bare Hamiltonian, the 



UV-divergent counter-term 



AW[r] = ~ / f\x-y\ 

Z Jx Jy 



-vqcL 



(D.7) 



which is a pure number, and thus does not change the expectation value of any physical 
observable. 

We next consider marginal operators: In the MOPE of Eq. (D.l), the integral over the 



relative distance of 



Jx—y \ 



-f- J -f- is logarithmically divergent at e — 0. In order to 

find the appropriate counter-term, we use dimensional regularization, i.e. set e > 0. An 
IR-cutoff L, or equivalently a subtraction momentum scale \i = has to be introduced 
in order to define the subtraction operation. As a general rule, let us integrate over all 
distances appearing in the MOPE-coeflicient, bounded by the subtraction scale L = /i -1 , 
giving 



+ 



+ 



L £ f(e,D) 



(D.8) 



Following Refs. [23-26,43], we use a minimal subtraction scheme (MS). The internal di- 
mension D of the membrane is kept fixed, and Eq. (D.8) is expanded as a Laurent series 
in e, which here starts at e^ 1 . Denoting by Res( | ), the residue of the term of order e^ 1 



of the Laurent expansion of ( 
to be 



for L — 1, the residue of the pole in Eq. (D.8) is found 



Res 



+ 



1 

2D 



(D.9) 



It is this pole that is subtracted in the MS-scheme by adding to the Hamiltonian a 
counter-term 

AH[r] = -^Res( 



(D.10) 



Similarly, the divergence arising from the contraction of two dipoles into a single dipole 
is subtracted by a counter-term proportional to the residue of the single pole of 



= / / 

J\x\<LJ\y\ 



•A<l 

J\x\<L J\y\<L V > 



■d/2 



(D.ll) 



Integrating the above yields 

Res/ 



L \2-D 



J -T(S) ' 

As a result, the model is UV-finite if we use the renormalized Hamiltonian 



(D.12) 



n R [r] 



2 - D Jx 2 



1 (Vr(x)) 2 + bZ bf i £ I I 5 d (r(x) - r(y)) , 



(D.13) 



xJy 
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instead of the bare Hamiltonian H[r]. Here, r and b are the renormalized field and coupling 
constant, and /i = L~ x is the renormalization momentum scale, and the renormalization 
factors at one loop order are 



Z = 1 - (2 - D)Res 
Zb = 1 + Res 



b 

e 



(D.14) 
(D.15) 



The renormalized field and coupling constants are re-expressed in terms of their bare 
counterparts through 



r {x) = Z 1 ' 2 r{x) , b, = bZ b Z d/2 ^ £ ■ 



(D.16) 



Following the analysis of Refs. [23-26] , the renormalization group /9-function and v (the 
anomalous scaling dimension of r) are obtained from the variation of the coupling constant 
and the field with respect to the renormalization scale /i, keeping the bare couplings fixed. 
They are written in terms of Z and Z b as 



1/(6) 



d_ 
dpi 



< HI 

bo 

2 — D 
2 

2 — D 



-eb 



l + bjL\nZ b + lb£ b \nZ 



(D.17) 



1 9 1 7 

-U-— InZ 

2 ^ „ 



(D.18) 



E Other topologies 

Throughout the manuscript we have generalized polymer loops at D — 1 to closed hyper- 
spheres at D 7^ 1. Natural questions are whether other topologies may also be used, or if 
it is possible to sum over all topologies with appropriate weights. The only equation where 
the topology enters, is that of the partition function for a single polymer or membrane in 
Eq. (5.11). Following David et al. [24], this is generalized to 



7 (o) 




+ - — 1 

D 6D 



(E.l) 



where x is the Euler characteristics of the membrane; x — corresponds to a sphere, x — 1 
to a 1-torus, and so on. Note that this equation is strictly correct only for D = 2, but 
that we make an analytic continuation for arbitrary D in order to keep the anomalous 
contribution from the topology. (Topological anomalies, related to trace anomalies in 
conformal field theory [70], occur only in integer dimensions.) Expanding now for e — 0, 
the 1-torus gives an additional contribution at D = 4/3, and d — 8. Thus for D > 4/3, 
the torus is irrelevant, higher topologies are even more irrelevant, and their neglect is 
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justified. In principle, for D < 4/3, we can perform a double e-expansion about this 
point. The second expansion parameter is 

6=e+±-D. (E.2) 

We can then introduce four different couplings. One coupling b on the same object (be 
it a torus or a sphere), a second coupling g between spheres, a third coupling u between 
spheres and tori, and a fourth coupling between tori. We leave such calculations for the 
future. 
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